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An extension of the Consistent-Q formalism for the Interacting Boson Model that includes the 
cubic {QxQx Q)^^-* term is proposed. The potential energy surface for the cubic quadrupole interac- 
tion is explicitly calculated within the coherent state formalism using the complete dependent) 
expression for the quadrupole operator. The Q-cubic term is found to depend on the asymmetry 
deformation parameter 7 as a linear combination of cos (87) and cos^ (87) terms, thereby allowing 
for triaxiality. The phase diagram of the model in the large N limit is explored, it is described the 
order of the phase transition surfaces that define the phase diagram, and moreover, the possible 
nuclear equilibrium shapes are established. It is found that, contrary to expectations, there is only 
a very tiny region of triaxiality in the model, and that the transition from prolate to oblate shapes 
is so fast that, in most cases, the onset of triaxiality might go unnoticed. 

PACS numbers: 21.60.Fw, 21. 10. Re, 05.30.Rt, 75.40.Cx 

I. INTRODUCTION 

The quadrupole operator plays a central role in Nuclear Physics because it is essential in the description of nuclear 
deformation, in the calculation of energy terms and in the evaluation of electromagnetic quadrupole transitions and 
moments [1 . Its presence is also of key importance within the Interacting Boson Model (IBM) |2l [3], where its 
components are defined, in terms of s and d bosons, as 

Q^ = [dt5 + ,tJ](2)+^[rftJ](2). (1) 

Various relevant operators (including scalar terms in the Hamiltonian or the E2 electromagnetic transition operator) 
are usually built from Q by appropriate tensor couplings. Within the IBM a useful and, at the same time, extremely 
simple Hamiltonian is the one of the Consistent-Q formalism (CQF) introduced some time ago by D.D. Warner and 
R.F. Casten [4 . This model Hamiltonian is formed by a one body monopole term, proportional to the number of 
d bosons, plus a quadrupole-quadrupole interaction among bosons. The CQF Hamiltonian can generate spherical, 
deformed axially symmetric as well as deformed 7— unstable potential energy surfaces [2 . It was noted [2^ that for 
a general s-d IBM Hamiltonian including up to two body terms, triaxial shapes are forbidden (please, note that 
triaxiallity can be obtained with up to two-body terms if g-bosons are included [5]). Within the standard s-d IBM 
several authors have shown (for example in Refs. [2|[6^^, in ch. 6 of Ref. [TO, ch. 2.2 of Ref. [11 and in Ref. ^12j ) 
that the triaxiality can be successfully introduced by adding three-body terms of the type [d^d^d^^^"^ ■ [ddd]^^\ These 
terms, together with a CQF Hamiltonian, generate a relatively broad region of triaxiality in the parameter space 
of the Hamiltonian. Although this is a valid way of generating triaxiality, it is not completely satisfactory because 
there is a priori no reason to invoke such terms and moreover one moves away of the desirable simplicity of the CQF 
framework. An alternative is to use higher order terms in Q as, for instance, the cubic {Q x Q x Q)^^^ interaction. 
The explicit expression for the cubic order interaction reads, 

(Q X Q X Q)(0) = ^(2^2 - ^|00)(2/i2(i/ - /i)|2^)g^Q,_^g_„ (2) 

where (....|..) stands for Clebsh-Gordan coefficients. Here three quadrupole operators (with rank-2 tensorial properties) 
are coupled to give a scalar term. This operator couples states with Ar = 3 as well as Ar = 1, where r is the 0(5) 
seniority quantum number. 

Several studies have been carried out for disentangling the properties of Hamiltonians containing Q x Q x Q terms. 
On the one hand, P. Van Isacker [13] has studied the tensorial properties of the Q x Q x Q operator, within the IBM, 
and has shown that, even with x = 0? this cubic operator can give spectra and band structures that are qualitatively 
similar to those usually associated with a SU(3) type of symmetry (axial rotor). On the other hand, Rowe and 
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Thiamova [14] have recently investigated the spectrum generated by using the fohowing Hamiltonian 

H = A^k{QxQx Q)(o) , (3) 

where A is the Casimir operator of 0(5), reaching similar conclusions to those of P. Van Isacker: the spectrum obtained 
by increasing the strength of the cubic term with x = displays the properties of an axially symmetric rotor. In their 
paper and, more recently in Ref. [15], it is shown that the contraction of the 0(6) algebra indicates that cos (87) is an 
image of the cubic term and therefore the only possible stable shapes are axially deformed. Following this reasoning, 
a quadratic term in the Q cubic scalar, i.e. {Q x Q x Q)*^^^ ■ {Q x Q x Q)'^^\ is necessary to generate potential energy 
terms displaying a cos^ (87) behavior. While this is most certainly true with their assumption on the form of the 

quadrupole operator (Q^ = [d^s + s^d]^^\ i.e. x — 0)? we find by explicit calculation that a term quadratic in the 
cosine is already present in the matrix elements of the {Q x Q x Q)^^^ operator alone, if one includes the x dependent 
term in the definition of the quadrupole operator as in Eq. ([T]). Therefore the possible onset of triaxiality can be 
studied at the level of the cubic CQF Hamiltonian (to be defined below) without the need to resort to higher order 
terms, in agreement with the already mentioned studies on the [d^d^d^^^"^ ■ [ddd]^^"^ operators. 

After this introductory section we introduce the cubic CQF Hamiltonian (Section [ll| and then obtain the potential 
energy surfaces (PES) generated by the {Q x Q x Q)^^'' operator. In Section [lll| we explore the phase diagram and 
discuss the order of the phase transition surfaces. Interesting limiting situations of this Hamiltonian are explored in 
Section IV and, finally. Section jVl is devoted to summarize our conclusions. 



II. CUBIC CONSISTENT-Q HAMILTONIAN 

The Consistent-Q formalism (CQF) is based on a simple IBM Hamiltonian [4 that allows to investigate not only 
the three dynamical symmetries of the IBM-1, but also the transitional regions in between. We propose in this section 
an extension of the CQF by adding to the original Hamiltonian a cubic combination of Q operators coupled to zero 
angular momentum, 

i^ = {..-(l-0[<%«+*,^'i4j«^l. (4) 

We call Eq. Q the Cubic Consistent-Q Hamiltonian (CCQH) because we consistently keep the same values of x all 
the quadrupole operators appearing in the above equation, up to the cubic degree. As usual, the various terms have 
been divided by the appropriate power of the boson number (N) in order to preserve the same kind of dependence 
in the large N limit and also to make sure that, in this limit, each term does not go to a infinity value. Although 
many other three body terms can be used (for instance the seventeen linear independent three body terms discussed 
in Refs. [H E] that have been successfully used in the interpretation of double phonon anharmonicities using the 
IBM), this looks like the simplest one and it is easy to justify on physical grounds as the first higher order interaction 
term in an expansion based on the quadrupole operator. The Q cubic term can be interpreted as a correction to the 
quadrupole-quadrupole scalar product rather than an additional term. 

This Hamiltonian has a rich structure that will be analyzed in depth in the following sections. Clearly, when ^ = 1, 
one falls back into the U(5) spherical limit. When ^ = and ks = one recovers the deformed 7— unstable (x = 0) 
and the axially deformed (x = ±a/7/2) IBM limits. For values ks ^ there will be a competition between different 
possible deformed shapes that could produce, in principle, stable triaxiality. In order to analyze this competition, and 
consequently the appearance of triaxiality, let us first investigate the geometry produced by the Q-cubic term in the 
Hamiltonian. 



A. Geometry of the {Q x Q x Q)^^^ operator 

One fundamental step in the investigation of the cubic Q term is the connection of the {Q xQ x Q)^^^ operator with 
geometry, that can be obtained within the well-known intrinsic state formalism. A coherent, or intrinsic, state [21 [3] 
is defined as a properly normalized application of the A^— th power of a linear combination of scalar and quadrupole 
boson creation operators to the vacuum, namely. 
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FIG. 1: (Color online) Left: surface (grid) and contours (colored) representing {Q x Q x Q)^'^'^ /N^ in the large N limit with 
X = as a function of {x,y}. The coordinates on the plane are x = /3cosj and y — /Ssinj. Right: the same function in the 
fundamental 0° < 7 < 60° wedge with contours (slightly different coloring scheme). 

where P and 7 are related to the deformation and asymmetry parameters, respectively. The main difficulty, when 
one deals with matrix elements of complicated operators within the coherent state formalism, is the length of the 
calculations. We have, therefore, set up a symbolic computer code that can evaluate in an analytic fashion complicated 
terms, keeping into account the non-commutative nature of the basic operators [16 . This code adopts the technique 
introduced in Ref. [17] of transforming the creation (annihilation) operators into derivatives acting on the left (right). 
After having carefully tested the procedure with several known quantities [10 , we have obtained the following result 
for the matrix element 

{/3,7,iV I (Q X Q X I /3,7,Ar) = , (6) 

i 

that is split for simplicity in one-, two-, and three-body terms (i = 1, 2, 3) given by 

To obtain the full matrix element these three terms have to be summed up, keeping in mind that, in the intrinsic 
state formalism, only the leading order in N is correctly evaluated and, hence, only ts is meaningful [T9]. 

From these expressions it is clear that the dependence on 7 when x = is of the type cos 37 and therefore the 
cubic operator with x = cannot generate triaxiality. The full expression, valid for any value of x, contains one term 
proportional to cos 37 together with another proportional to cos^ 37. Therefore these contributions, in addition to a 
CQF energy surface, might generate a triaxial minimum in the potential energy surface. 

In order to analyze the energy surface (|6|, we will examine in Fig. [l] the potential energy surface given by 
{Q X Q X Q)'^^^/N^ in the large N limit for x = 0. The upper left-part of the figure displays the value of the 
matrix element on the vertical coordinate, while the lower left-part is its projection on the plane with coordinates 
X = /3cos7 and y = Psmj. The surface has prolate minima at /3 = 1 and 7 = and oblate maxima at /3 = 1 and 
7 = f (1 + 2n) with n = 0, 1,2. The contour map on the x — y plane gives a color-coded projection of the surface 
values, with black/blue representing the (prolate) minimum and yellow/red representing the (oblate) maximum (see 
color bar on the right). For completeness sake we also give in the right-part of the figure a contour plot in the more 
standard (/3, 7) polar coordinates, limited to the 0° < 7 < 60° wedge with a slightly different coloring scheme. Of 
course the character of the minima/maxima can be interchanged by an overall sign change. In conclusion, the cubic-Q 
term with x = o^ly produces either prolate or oblate minima. Triaxiality or 7— unstability are not allowed with the 
Q X Q X Q term if x = 0. 
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III. THE PHASE DIAGRAM IN THE LARGE N LIMIT 

The study of the shape of a quantum system, in particular of atomic nuclei, proceeds through the study of the 
energy surface in the large N limit. The shape is strictly defined in the thermodynamical limit. The procedure starts 
by calculating the potential energy surface in the large N limit, then a minimization should be done, getting the 
equilibrium value of the deformation parameters. 



A. Potential Energy Surface 

The energy per boson in the large N limit can be readily calculated from Eq. (|4|): 



J?(C,fc3,X,/?,7) = C 
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(10) 



Note that the N dependence is eliminated due to the proper N scaling of the different terms in the Hamiltonian and 
to the use of the energy per boson. 



The presence of cos 87 and cos^ 87 in (10) can give rise to triaxiality. Let us start discussing the possibility of 



triaxiality for a simpler situation: an energy surface independent of /3, but with the 7 dependence as in Eq. (10). In 
that case, one can write the energy surface as 

£^(7) = acos^ 87 + 6COS87 + c (11) 

where we take a, 6, and c as constants for the moment. This function admits triaxial extrema at 

7 = (± arccos {-b/2a) + 2n7r) /8 (12) 

if I —b/2a\ < 1, with n = 0,1,2. When instead | —b/2a \> 1, the extrema sit at 7 = 27rn/8 and 7 = (2n + l)7r/8, with 
n = 0, 1, 2, and triaxiality is excluded. 

In the general case, in which the coefficients a and b depend on /3, as in the energy functional Eq. ([lO|, the 
presence of terms depending on cos 87 and cos^ 87 can, in principle, produce triaxial shapes. It is clear that the 
existence of a triaxial minimum is linked to the inclusion of the cubic-Q term {ks 7^ 0) in the IBM Hamiltonian. The 
argument discussed above could, in principle, be valid even when the coefficients a and b depend on /3, as long as 
their dependence on the equilibrium value of /3, /3o, is smooth enough with respect to the control parameters of the 
Hamiltonian. However, for the cubic CQF Hamiltonian the latter argument is not so obvious, mainly due to the non 
continuous or very "fast" dependence of the P equilibrium value in some intervals of the control parameters. To shed 
some light on this, one can calculate the value of | — 6/2a | (see Eq. ([T2|)), which is the only necessary quantity to be 
known for getting the value of 7, and it can be written as 



b 




~2a 





(13) 



where f3o is the equilibrium value of /3— , that is, at the minimum of the energy. To rule out the possible existence of a 
triaxial region one has to prove as necessary and sufficient condition that | ^ | > 1 , for any value of the parameters in 
the Hamiltonian. In fact, this is true for ^ = 0. Note that Eq. (18) does not depend explicitly on ^, but its dependence 
comes in implicitly through the dependence of the equilibrium value of /3 on ^, x, 





and ks. 

To start with, we can easily calculate the value of |^| in some limiting situations. In particular, lim^^oo \ ^\ 
or lim|^|^oo \ ^\ — ^ therefore 7 = 7r/6, but these situations are quite unrealistic. More interesting limiting 
situations are: lim/3^0 \ ^\ = co, lim/cg^oo \ ^\ > 16 (this case is obtained considering that fS^x^ is always larger 
than zero, remembering that /3o in this equation is the value at the equilibrium point, taking into account that the 
larger is |x|, the larger is f3o and, finally, fixing the maximum values for |x| and f3o to a/7/2 and l/\/2, respectively). 



limz, 



00, and lim^ 



00. We will also prove later on that 
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> 1 for ^ = 0. On the other hand 



it is also true that the value of /3 obtained for ^ = is always larger than the corresponding one for ^ 7^ (for the 
same values of x and ks). With these general ideas in mind, it is easy to see that |^| >> 1 except, maybe, for a 
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very narrow range of the parameters. Indeed, one can see that, once the values of the Hamiltonian parameters have 
been set, if we treat ^5 as a free parameter and not as a equihbirum value, one gets |-^| < 1 only in a very narrow 
range of /3. In particular, it is possible to see, numerically, that this particular range of P is placed in a region where 
/3 exhibits a sudden drop. This is illustrated in Fig. [2j where the value of | ^ | is plotted as a function of /3, for /ca = 1 
and X = —0.5. It is clearly observed that only around P ~ 0.934 one gets |^| < 1, but this particular value of P falls 
in a rapidly changing region as it can be appreciated in the lower part of the figure, where the equilibrium value of {3 
is represented as a function of ^. 




I I [ I I I I 

0.2 0.4 0.6 



FIG. 2: Upper part: value of | (see text) as a function of (3 for k^, — 1 and x — —0.5. Equilibrium value of /3 as a function 
of^. 



As a conclusion of this qualitative discussion on possible triaxiality produced by the Hamiltonian Q, it can be said 
that one expects that the region supporting triaxial shapes is small, if any. We will quantify this idea in the following 
sections. 



B. The coordinates 



The pictorial representation of the phase space of the CCQH is an extension of the well-known Casten triangle 
(left plot in Fig. [3|. In our case the parameter space becomes a tetrahedron, where the horizontal coordinates are 
related to the ^ and x control parameters in the Hamiltonian, while the vertical coordinate is directly connected to 
the coefficient of the Q x Q x Q term. We use the following parametrization: 

p=l-i; ^=-'^ z = pks]. (14) 



3 V7 

In the studies presented below, the value of x ranges from — >/7/2 to and therefore the amplitude of the (j) angle is 
30°, p ranges from to 1 and ks is taken as positive. This range of control parameters allows to get the results for the 
full model space since the energy surface (10) presents the following prolate/oblate symmetry: x ~^ ~Xi ^3 ~^ ~^?>i 



7 ^ 7r/3 — 7. This allows to extend the results obtained in the case of x < and /C3 > to the regions with x > 
and k^ < 0. 

It is worth noting that due to the form of the Hamiltonian Q it is not possible to reach a pure Cubic Q Hamiltonian, 
except for ks ^ 00. However in Fig. [3]the upper part of the tetrahedron is labeled with Q xQ x Q to denote that this 
limit is reached for large values of z = pks. Anyway, in this study we will limit ourselves to moderated large values 
of ks (up to 10) since the cubic Q term is supposed to be a correction to the dominant QQ term. 

First, the schematic phase diagram of the model Hamiltonian will be studied and then some particular relevant 
regions will be analyzed in more detail with the aim of singling out possible regions of triaxiality. 
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FIG. 3: (Left) Schematic representation of the CCQH parameter space and coordinates. Four quadrants are present, of which 
only the first one is studied in detail as explained in the text. As a reference, the IBM dynamical symmetry limits are explicitly 
shown. (Right) Schematic phase diagram for the CCQH, corresponding to quadrant I in the figure on the left. 



C. The geometry of the phase diagram 

The resulting phase diagram for the Hamiltonian (|4) is depicted in the right hand side of Fig.js] Although a detailed 
explanation of the nature of this phase diagram will be given along this section, we already present the complete 
diagram here to facilitate the following analysis. The position of the critical surfaces are determined numerically, 
except for some particular regions of the phase diagram that can be obtained analytically. We have explored the 
parameter space to find out the position of the critical surfaces, following certain selected paths across the parameter 
space to clearly illustrate where the critical surfaces are placed and what is their character. This phase diagram is 
characterized by the existence of four different phases: spherical (p < 0.2, or 1 > ^ > 0.8), prolate and oblate axially 
deformed shapes {p > 0.2, or ^ < 0.8) and a very tiny region of triaxial shapes, occurring between the oblate/prolate 
separation surface close to the x = ±>/7/2 faces of the tetrahedron. Note that the existence of a spherical region is 
limited to moderate values of ks. Due to the particular form of Hamiltonian Q, very lare values of ks will transform 
the spherical region almost in a single point around the origin. The two plotted surfaces correspond to first order 
phase transitions except along their intersection line, which is a second order phase transition line (full line in the 
right panel of Fig. [3|. As we will explain in detail along this section, the prolate-oblate first order transition surface 
becomes a two-fold second order phase transitions surface for x ^ — although these surfaces are extremely close 

and cannot be distinguished in Fig. [3] 

The full phase space comprehend four tetrahedra (numbered I - IV) with a common U(5) vertex, shown in the 
left part of Fig. [sj while in the right part of the figure only the upper left quarter (I) with positive values of ks and 
negative values of x is plotted. Apart from the spherical region close to ^ = 1, that is always present, the upper right 
quadrant II {ks > and x > 0) contains only oblate shapes, while the opposite lower left quadrant HI {ks < and 
X < 0) contains only prolate shapes. Finally, due to the prolate oblate/symmetry x ^ — ^3 ~^ —ks^ 7 ^ 7r/3 — 7, 
the lower right quadrant IV can be obtained by rotating the first quadrant 180° with respect to the U{5) — 0(6) axis. 
In order to illustrate how the phase diagram has been obtained we will show the results along several selected paths 
within quadrant I. 



1. Paths from deformed to spherical shapes 



First, we will explore the separation surface of spherical and deformed shapes. For that purpose, we start selecting 
/cs = 4 and x = —a/7/2 and we vary ^. In Fig. we observe a first order phase transition, where /3 goes from a finite 
value to zero and 7 goes from 60° (oblate shape) to undefined (spherical), when entering the spherical region. Note 
that the equilibrium value of /3 for ^ = is /3 = 1/ >/2. In the right part of the figure we have displayed the position of 
the path (green thick line) inside the phase diagram and also the position of the phase transition point (green circle) . 
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FIG. 4: (Color online) Ground state energy and the equilibrium value of the shape variables /3 and 7 as a function of ^ along 
the path shown in the right part of the figure, for /cs = 4 and x — ~V7 /2. Notice that the equilibrium value of 7 is undefined 
in the spherical region. 




FIG. 5: (Color online) The same as in Fig. [i] but for = 1 and x — — a/7/2. 

It should be noted that U(5) corresponds to ^ = 1 and the surface SU(3)-0(6)-QQQ to ^ = 0. In the following figures 
we include plots with similar interpretations. For x = one can obtain plots similar to Fig. |4j although in this case 
the transition occurs for a slightly larger value of ^. Now we consider a smaller value of k^. We impose ks = I and 
X = —V7/2 and we vary the value of ^. This trajectory is plotted in Fig. [H] and once more the existence of a first 
order phase transition can be noted, although it is not as abrupt as in the previous case. The value of 7 passes from 
zero (prolate shape) to undefined in the spherical region. Finally note that /3 = V2 for ^ = 0. 

Now, it is clear that there is a surface that separates spherical and deformed shapes at around ^ = 0.8. In order to 
investigate the character of this surface, where the system changes from a spherical to a deformed shape, it is possible 
to carry out a Taylor expansion of the energy around the value /3 = 0, obtaining: 

E{^,ks^X.Pn) = (-4 + 50/3' 

- ^ (1 - (2 VTOfe + 5 vTIx) cos 37 

+ (-e + (l-0(8-^-'f))/3^ + 6(n (15) 



FIG. 6: (Color online) The same as in Fig.raor|^but for = 1 and x — —2/ a/5 (second order phase transition). 



This expression is very convenient because one can easily read off the order of the phase transition when crossing the 
surface. In general, the presence of a cubic term in j3 implies that the system undergoes a first order phase transition 
[20] , while its absence guaranties that the phase transition is of second order. 



From equation ( 15 ) and the preceding discussion one notices that the spherical-deformed surface corresponds to a 



first order phase transition except when 



fc3 = -^X, (16) 



which cancels the term and changes the transition type to second order. This second order phase transition line is 
the intersection of the spherical-deformed surface (studied here) and the prolate-oblate surface (to be studied in the 
next subsection). Therefore along this line spherical, prolate and oblate shapes coexist. 

In order to show the characteristics of this second order phase transition line we have followed a path that crosses 
the spherical-deformed surface precisely through this line. This is illustrated in Fig. [6| where a second order phase 
transition appears at ^ = 4/5. No discontinuity is observed neither in the order parameter /3, nor in the energy or 
in its first derivative (not shown in the figure). The discontinuity appears in the second derivative of the energy 
confirming that the line (p^Gl) is of the second order type. 



With respect to the existence or not of triaxiality in the region close to spherical shapes, in expression (15) only 
the cos 37 shows up, but the cos^ 87 dependence will appear together with {3^ terms and therefore it will generate a 
value >> 1 since the coefficient of cos 87 (a) is much smaller than the one of cos^ 87 (b) for /3 < 1. Therefore, 
there is no possibility of triaxiality close to the spherical region. 

Finally, note that for every value of ^, with the exception of ^ = 1, it always exists a value of that marks the 
appearance of an additional deformed minimum. This situation is quite different with respect to the one obtained 
with two body Hamiltonians, where a value of ^ > 0.8 implies a spherical minimum regardless of the value of x- 
the phase diagram depicted in the right part of Fig. [3] this fact cannot be noticed because of the moderated values of 
/c3 used along the diagram. 



2. Paths from deformed prolate to deformed oblate shapes 

Now, we explore the separation surface between prolate and oblate shapes, starting from the ^ = surface and 
then generalizing to^^O. 

a. The ^ = case: interplay of Q ■ Q and {Q x Q x Q)^^^ 

In this subsection the ^ = (p = 1) surface (SU(8)-0(6)-QQQ) is studied. For pedagogical reasons we 
start our study by setting x = —V7/2 (the vertical line between SU(8) and QQQ in Fig. s]) and later on we will 
generalize the outcome to arbitrary values of x- For the case, ^ = 0, x = —a/7/2, the extrema of the —Q • Q/N'^ term 
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FIG. 7: (Color online) Ground state energy, equilibrium shape variables 13 and 7 as a function of fcs for ^ = and x = — \/7/2. 
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d'^E/dp'^ d'^Ejd-i^ character 







- M 




7r/3 


+ - sp 







+ + m 


2V2 


7r/3 


- - M 



TABLE I: Extrema of -Q • (5/A^^ for x = -^7/2, along with their character. The sign of the second derivatives with respect 
to /3 and 7 are also given (see text). 

in the large N limit are listed in Table |lj In this Table the signs of the second derivative of the energy with respect 
to both shape parameters, /3 and 7, and the character of the extrema (M for maximum, m for minimum and s.p. for 
saddle point) are given. To avoid equivalent shapes related with the symmetries of the P.E.S. 

we impose the constraints /3 > and < 7 < 7r/3. Summarizing the information in Table |lj the surface obtained 
from the —Q • QjN^ term presents: a spherical maximum with no dependence on 7, a prolate minimum for /3 = V2, a 
saddle point in the oblate side for /3 = 1/a/2 and an oblate maximum for /3 = 2^/2. Note that at the extrema there are 
no off-diagonal terms in the Hessian matrix. A similar study can be done for — (Q x Q x Q)^^^ jN'^ . The characters of 
the extrema for the surface produced by this term are listed in Table [TH In this case there are a spherical maximum, 
an oblate minimum at /3 = 1/a/2, a prolate maximum at /3 = a/2 and, finally, an oblate saddle point at /3 = 2\/2. 
Moreover there is also a minimum for /3 ^ 00 and 7 = 7r/6. Note that, again, at the extrema there are no off-diagonal 
terms in the Hessian matrix. It is worth noting that the extrema for both, —Q ■ Q/N'^ and —{Q x Q x Q)^^'^/N^ 
coincide, although their character, either minimum or maximum, can be different. 

A Hamiltonian built upon a linear combination of —Q • Q and —{QxQx Q)*^^^ is expected to possess eigenstates 
with large quadrupole moments and small fluctuations around their equilibrium values. This kind of state is eigenstate 
of the quadrupole operator, as well as of both, Q ■ Q and {Q x Q x Q)^^^ [3l[T2]. This is the reason why the extrema 
points of both parts of the Hamiltonian under study coincide [12j. 
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s.p. 



TABLE II: Same as Table [l] but for -(Q x Q x Q^^^/N^. 
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FIG. 8: (Color online) Energy surface for -Q • Q/N'^ - x Q x Q)'^^^ /N^ with x = -^7/2, corresponding to the green 

dot of the right part of the figure. Note the existence of a prolate and an oblate degenerated minima. 
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TABLE III: Position of the possible minima for a linear combination —Q • Q/N'^ — k3{Q x Q x Q)^'^'^ /N^, with a generic value 
of X < 0. 

In general, a linear combination of —Q • Q and —{QxQx Q)*^^^ always generates a spherical maximum for (3 = 0, 
an oblate maximum or saddle point at /3 = 2^/2 and a competition between two minima, one prolate at /3 = a/2 and 
the other one oblate at /3 = 1/ V2. The above situation corresponds to a first order phase transition and it excludes 
the presence of a triaxial minimum. 

Fig. [7| gives the evolution of the ground state energy and the equilibrium value of f3 and 7 for the path along the 
line X = —a/7/2 and ^ = 0. This figure shows the typical ingredients of a first order phase transition. There is a 
discontinuity in the first derivative of the energy and a discontinuity in the value of f3 and 7. For low values of ks 
the minimum corresponds to f3 = \/2 and 7 = 0°, while for large values to /3 = 1/a/2 and 7 = 60°. The critical value 
of the control parameter corresponds to kz = >/35/3. This figure is similar to Fig. 2 appearing in Ref. [21] although 
there the phase transition is second order and 7 goes from 0° to 60° in a smooth way and therefore a broad region of 
triaxiality exists. 

An interesting case happens for k^ = for which the prolate minimum (/3 = a/2) flattens (becoming a maximum 
for larger values), i.e. d'^E/dp'^ = d'^E/dj'^ = 0. Additionally, the flat prolate and the oblate minima become 
degenerate with energy £^ = — |, being also degenerate with the minimum at ^5 ^ 00, 7 = 7r/6. This is illustrated in 

Fig. [s] For ks < the only minimum is the prolate one, that becomes flat for ks = For ks > the only 

minimum is the oblate one. Therefore, at ks = the spinodal and the antispinodal points coincide and also are 
degenerate (first order phase transition point). 

It is possible to extend the preceding study to a generic value of x- Without loss of generality, we choose X < 0- 
this case the position of the possible minima are given in Table |III[ Again there is a competition between the prolate 
minimum, coming from the —Q ■ Q/N'^ term an the oblate one, generated with the —{Q x Q x Q)^^'^/N^ term. In 
this case there exists a coexistence region, but again there is no room for triaxiality. In particular, the position of 
spinodal (7 = 7r/3) and antispinodal (7 = 0) points are given by 
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49\/5x (x' ± xVx^ + 14 + 14 



± 4x5 + 14 + eSx"^ ± 35x3 + 14 + MTx^ + 686 



(17) 



the upper (lower) sign corresponds to 7 = 7r/3 (0). The value of ks where the two minima become degenerate (critical 
line) reads 



he = 



2x^ + 7 



(18) 



In Fig. [9] the spinodal, the antispinodal and the critical values of ks are represented in the {ks^x) plane for ^ = 0. 
This figure clearly shows the existence of a coexistence region and the absence of a triaxial region in the ^ = surface. 
Additional extrema to the ones given in the above mentioned table can appear which do not correspond to minima. 




FIG. 9: (Color online) Value of the spinodal (7 = vr/S), antispinodal (7 = 0) and critical value lines of as a function of x for 



b. The ^ 7^ case: The prolate- oblate-triaxial critical surface 

The analytical arguments used in the preceding discussion are no longer valid when the U{5) term gives a 
contribution to the energy surface, i.e. 0, however they suggest that the parameter space region that gives rise to 
a triaxial minimum should be rather small. It is not possible to obtain analytically the expression of the equilibrium 
value of 7 for a general energy surface. 

To study numerically the deformed region for this particular Hamiltonian we calculate the value of the energy, /3 
and 7 along some particular lines that go through the phase diagram and cross the prolate-oblate separation surface. 
In Fig.[lO]we select a particular path with ^ = 0.5 and x = — 1 with kz ranging from to 4. In this figure the presence 
of a first order phase transition shows up because there exists a clear discontinuity in the values of /3 and 7, when 
crossing the prolate-oblate surface. Therefore there is a jump from the prolate to the oblate minimum at the point 
where they become degenerate and the triaxiality is forbidden. The same happens for other values of ^ and x- 

However a more detailed inspection shows that there exists a very tiny region of triaxiality around x ^ ~V^/2, for 
^ < 4/5. This is clearly illustrated in Fig. 11, where the lines that separate the prolate, oblate and triaxial regions 
are plotted for a value of ^ = 0.5. Note that the triaxial region occurs for x close to the limiting value of — >/7/2. 
At this limit k^ G [1.752 — 1.791]. The lower part of the figure corresponds to prolate minima, while the upper part 
to oblate minima. The box on the vertical left side shows the blowup of this region: there is indeed a small range of 
values where it is possible to find triaxial minima by numerical procedure. The more one moves away from the value 
the narrower this range becomes, due to the delicate interplay between the quadratic and cubic terms. 



X 



Here we have plotted several potential energy 
with X — ~^ and ^ = 0.5. The minimum is 



We can see in a clearer way the onset of a triaxial minimum in Fig.^ 
surfaces corresponding to points along the left vertical axis of Fig. Ill 
clearly prolate for small values of /C3 and it starts to widen as one approaches the region k-^ ^ 1.7. From approximately 
1.752 to 1.791 the minimum is triaxial, but very shallow and therefore invisible at the present scale of 7). We show 
in Fig. [13) for completeness sake, the close up of the minimum for /C3 = 1.77. Increasing again the value of /C3 leads 
to an axially deformed oblate minimum. It is worth noting that the change through the triaxial region is extremely 
swift and some care must be taken in the minimization. One can appreciate the extremely shallow character of the 
triaxial minimum that extends from the axially prolate till the oblate minimum, very much in a similar way to Fig. [S] 
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FIG. 10: (Color online) Ground state energy and equilibrium value of the shape variables /3 and 7 as a function of fe, for 
^ = 0.5 and x — ~^ for the path shown on the right. 




FIG. 11: (Color online) Section of the three-dimensional parameter space with ^ = 0.5 and N 00. The yellow inset shows 
the tiny triaxial region (see text) that is found at x ~ —a/7/2. 



where prolate and oblate minima are degenerate and no triaxial minimum exists. The existence of the triaxial region 
is confirmed by inspecting contour plots of the potential energy surface along a suitable trajectory of the parameter 
space, rather than just looking at the plots of /3 and 7 as a function of the parameters ^, x ^3- 



Finally we present in Fig. 
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a cut of the triaxial region as a function of ^ and ks for a constant value of x = —V^ /2. 
A very small triaxial region is always found, except for ^ = 0, between the lower prolate region and the upper oblate 
one (see inset of Fig. 14) and terminates in correspondence with the spherical phase that occurs for high values of 
^. In a preceding subsection we have discussed that the point connecting with the spherical phase is indeed a single 
point as well as the one at ^ = 0. 

A final question should be answered: what is the order of the phase transition developed in crossing the prolate- 
triaxial-oblate surface? In the region of triaxiality the surface evolves into two different sheets where second order 
phase transitions exist. As far as we move from x = —V7/2 this two-fold surface collapses in a single one and the 
phase transition becomes of the first order. 



Since the small triaxial region appears in the deformed area close to x = we present in Fig. 15 a path with 
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60° 



50° 



50° 




k3 = 0.0 



k3 = 1.5 



k3 = 1.77 



k3 = 2.0 



k3 = 3.0 



FIG. 12: (Color online) Potential energy surfaces for the CCQH for x — —a/7/2 and ^ = 0.5, with different values of 
indicated in the figure. The observed minima range from an axially deformed prolate minimum for /cs = 0.0 and = 1.5 to 
an oblate one for = 2.0 and = 3.0 , passing through the triaxial region. The triaxial minimum when = 1.77 is not 
apparent from this figure, because it is very shallow (see next figure). 



60° 




FIG. 13: (Color online) Closeup of the potential energy surface of Eq. (10) with x — —a/7/2, <J = 0.5 and = 1.77, clearly 
showing a shallow triaxial minimum at /S ^ 0.62 and 7 ^ 23.4°. One can prove, with similar plots, that the prolate minimum 
smoothly shifts through the triaxial region and finally becomes oblate. 



fixed X = — V^/2 and ^ = 0.5 and changing k^. Once more we plot in this figure the evolution of the ground state 
energy and the equilibrium value of /3 and 7. Again low values of ks give rise to prolate shapes, while the larger ones 
produce oblate forms, but a narrow region of triaxiality exists around ks ^ 1.75 — 1.79. The scale does not allow to 
discriminate, therefore we have enlarged it in the inset that shows two, very close, second order phase transitions. 
One from prolate to triaxial at around = 1.75 and a second one from triaxial to oblate at around /cs = 1.79. As 
soon as we depart from the x = —V7/2 surface these two second order surfaces approach more and more up to a 
point in which both coincide and apparently transform into a first order surface that directly separates prolate and 
oblate shapes, without any triaxiality in between. It is very difficult numerically to determine whether the triaxial 
region narrows indefinitely or rather ends up in a tricritical point. 
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FIG. 14: (Color online) Section of the three-dimensional parameter space with x — —^ as a function of ^ and fe. The inset 
shows the triaxial region in yellow on a finer vertical scale (see text). 




FIG. 15: (Color online) Ground state energy, shape variables ^ and 7 as a function of A:3 for ^ = 0.5 and x — —a/7/2 for the 
path shown on the right. 



3. A path going from spherical through oblate and prolate shapes 



As a final calculation we study a trajectory that crosses the spherical-deformed surface as well as the prolate-oblate 
one. The appropriate values of ks and x constrained by equations (IT6| and (18) and should verify 



2x^ + 7 



> h > 



V5 



-X- 



(19) 
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and 



We choose ks = 1.5 and x = — 2/>/5, which fulfil ( |19| ), and vary the value of ^. This path is depicted in Fig. 
two first order phase transitions are clearly marked. From left to right, the first one has a discontinuity in (3 between 
two finite values changing at the same time 7 from 0° to 60°. In the second one the system goes from the oblate 
phase to the spherical one, changing (3 from a finite to a zero value and 7 from 60° to undefined. 
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FIG. 16: (Color online) The same as in Fig.jijbut for = 1.5 and x — — 2/a/5. 

IV. SOME APPLICATIONS 
A. Spherical to axially deformed critical point: X(5) 




-60° 



FIG. 17: (Color online) Potential energy surface (left) obtained from (|4| with x = 0,^ — 0.9 and = —9.5. The valley 
along the 7 = 0° direction for prolate deformations resembles the potential used in the X(5) solution of the Bohr Hamiltonian 
[22]. Right hand side: two sections of the energy surface around the equilibrium value of the deformation parameters have been 
plotted as a function of (3 and 7, respectively. 



An interesting special case is to keep x — ^-> order to eliminate the driving effect of the Q • Q term toward axial 
deformation. In the left part of Fig. [T7|a large value of k^, generates a surface with an almost flat valley along the P 
direction, while the prolate (oblate) confinement in 7 is due to the negative (positive) sign in front of the {QxQx Q)*^^^ 
term. This case has a potential energy surface that has a flat shape qualitatively similar to the V(/3,7) potential 
used in the X(5) critical point symmetry [22], with due differences: it is generated through a set of parameters X = 0? 
^ = 0.9, and ks = —9.5 (a positive value would generate an equivalent oblate shape), the potential does not tend to 
infinity, but goes smoothly to an asymptotic value and the periodicity in 7 is retained, at a variance with a harmonic 
oscillator. The bottom of the potential is not exactly flat, nor the behavior approximates a dependence, but the 
fact that a potential energy surface that mimics X(5) could be obtained with x = is somewhat surprising, because 
it has been used to associate X(5) to a case that is intermediate between spherical and axially prolate shapes, while 
here we do not even need to set x = —V7 /2. 
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60° 60° 60° 60° 




k=0.0 k-0. 58776 k=1.0 k=10.0 

FIG. 18: (Color online) Potential energy surfaces for the Hamiltonian ([3| with = 10 and different values of k : 0.0, /ccr, 1-0 
and k — 10.. The minimum (blue) ranges from spherical to axially deformed prolate. Vertical scales change along the transition. 



B. Application to a schematic Hamiltonian 



As an application, we can now use the formula for the matrix element of the (Q x Q x Q)^^^ operator with X = 
within the boson coherent state formalism to give further insight into the results obtained by Rowe and Thiamova 
for the schematic Hamiltonian ([3| [HI [15] . The parameter k is allowed to vary from zero to large positive values. 
In Fig. [18] we have plotted several contour maps corresponding io k = 0.0, 0.58776, 1.0, and 10.0. Notice that the 
A term alone {k = 0) just gives a spherical minimum, because its matrix element in the coherent state approach 
is just proportional to From the figure one can appreciate how, with growing value of /c, the minimum moves 
from the spherical configuration to the axially deformed prolate one. This is due to the competition between the 
spherically driving term A and the cubic interaction term (Q x Q x Q)^^^ that has a stable axial minimum, despite 



having x = 0. We have collected in Fig. 19 several cuts of the potential energy surface along 7 = 0°, for the values 



k = 0.0, 0.2, 0.4, kcri 0.8, and 1.0. Alongside the ever-present minimum at /3 = 0, this potential admits a second 
deformed local minimum, that starts to appear after k ^ 0.51. The critical point for the first order phase transition is 
found around kcr ~ 0.58776, where two minima are degenerate. After this point the deformed minimum prevails. The 
values discussed above are valid with the particular choice of A" = 10 and with the Hamiltonian as given originally 
in ([3|. Although it is immaterial to the present discussion, one should notice that, while the 0(5) scalar term is 
linear with the boson number, the {Q x Q x Q)^^^ term is cubic with N and therefore one should better choose a 
Hamiltonian that is properly normalized to get rid of the different scaling with the boson number. Notice that Ref. [13] 
uses a slightly different Hamiltonian having L • L as a spherical term: although it is most certainly important to use 
the proper terms when calculating an IBM spectrum, it is less critical here, since both these operators have matrix 
elements proportional to (albeit not with the same coefficients). In both cases [131 [Mj triaxial minima are clearly 
not possible. Moreover, the inclusion of the cubic term, even with X = 0? forbids the existence of a 7— independent 
solution. 



V. CONCLUSIONS 



We have calculated the potential energy surface of the cubic {QxQx Q)^^^ term within the coherent state formalism 
with the most general expression for the quadrupole operator. This has allowed us to confirm the results of Ref. [13] [T4] 
concerning the fact that this term with x = can generate axially deformed minima, but in addition it has been 
shown that the cubic term has already a dependence on cos^ 37 (for X 7^ 0) together with the dependence on cos 37 
(obtainable also when X = 0)- Therefore a Hamiltonian containing this term may generate triaxiality without the 
need to resort to the more complicated {Q x Q x Q)'^^^ ■ {Q x Q x Q)*^^^ expression of Ref. [T5] . 

After a discussion of its general features, we have used it to extend the Consistent- Q Hamiltonian to the Cubic 
Consistent-Q Hamiltonian, of which we have discussed the phase space, discovering that a triaxial region can indeed be 
found, between the oblate and prolate phases. This region is extremely tiny (at least in the particular parameterization 
chosen, if compared to the other phases), and can be pinpointed numerically only very close to the limiting values of 
X ^ ± a/7/2. The inspection of the PES confirmed that the minimum is indeed triaxial inside this region. It is not 
clear, at the moment, if the prolate and oblate phases are always separated by a region of triaxiality that progressively 
tails off as one goes to x = or if the triaxial region disappear at some point. The numerical results suggest that this 
tiny triaxial region finishes in a line where prolate, oblate and triaxial shapes coexist, i.e. a tricritical line. 
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FIG. 19: (Color online) Cuts along 7 = 0° of the potential energy surfaces for the Hamiltonian ^ with x — ^ = 10 

for several values of k from 0.0 to 1.0. At the value kcr ^ 0.58776 it shows the typical signs of a first order phase transition. 
Purple dots indicate the position of minima. 
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